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ABSTRACT 

Context. It has been proposed that mergers induce starbursts and lead to important morphological changes in galaxies. Most studies 
so far have focused on large galaxies, but dwarfs might also experience such events, since the halo mass function is scale-free in the 
concordance cosmological model. Notably, because of their low mass, most of their interactions will be with dark satellites. 

Aims. In this paper we follow the evolution of gas-rich disky dwarf galaxies as they experience a minor merger with a dark satellite. 
We aim to characterize the effects of such an interaction on the dwarf’s star formation, morphology, and kinematical properties. 
Methods. We performed a suite of carefully set-up hydrodynamical simulations of dwarf galaxies that include dark matter, gas, and 
stars merging with a satellite consisting solely of dark matter. For the host system we vary the gas fraction, disk size and thickness, 
halo mass, and concentration, while we explore different masses, concentrations, and orbits for the satellite. 

Results. We find that the interactions cause strong starbursts of both short and long duration in the dwarfs. Their star formation rates 
increase by factors of a few to 10 or more. They are strongest for systems with extended gas disks and high gas fractions merging 
with a high-concentration satellite on a planar, radial orbit. In contrast to analogous simulations of Milky Way-mass galaxies, many 
of the systems experience strong morphological changes and become spheroidal even in the presence of significant amounts of gas. 
Conclusions. The simulated systems compare remarkably well with the observational properties of a large selection of irregular dwarf 
galaxies and blue compact dwarfs. This implies that mergers with dark satellites might well be happening but not be fully evident, 
and may thus play a role in the diversity of the dwarf galaxy population. 

Key words. Galaxies: dwarf - Galaxies: evolution - Galaxies: interactions - Galaxies: irregular - Galaxies: starburst - (Cosmology:) 
dark matter 


1. Introduction 

In the lambda cold dark matter (ACDM) paradigm, small dark 
matter halos are abundant. Most of these halos, with Myir < 
IO^Mq, are predicted to be strongly affected b y reionization . 


photo-evaporation, and/or supernova feedback (iGnedin 


Hoeft et ^ 2006t Kaufmann et al.l I2007 : Okamoto,et^ 


200C; 


2008 


Gnedin et al.l 120091: iLi et al.H2010l: Sawala et al. j2013l: but see 
also iTavlor & Webstej 2005t IWarren et al. 20071) . These pro¬ 
cesses thus cause progressively larger numbers of small dark 
matter halos to be almost completely dark. The existence of such 
dark galaxies is a solution that is often suggested to the m issing 
satellites problem dKlvnin et al.lll999l:lMoore et alJll999l) . 

Dwarf galaxies are known t o be very inefficient a t forni - 
ing stars (iBlanton et alJ 120011 : [Robertson & Kravtsov 200^ 
to have very low baryo n fractions ( GnedinI l2000l: Hoeft et alJ 
I 2 OO 6 I: ICrain et ^l200^ . and to generally be gas-rich if in the 
field. This is consistent with the expectation that the stellar- 
to-halo-mass ratio must decrease steeply toward lower masses 


( Behrooz i et al.l20'T3: Moster et al.l2013l:lKormendy & FreemanI 


must have much lower baryon ic component mass es, and so most 
of their satellites will be dark dHelmi et al.ll2012h . 

Mergers of gas-rich galaxie s are often thought to gi ve 
rise to bursts of star formation dMihos & Hernaui^ Il994alf5 : 
iTevssier et al.l 1201(1 iBournaud et al.l 1201 ll) . although simula- 
tions suggest this de pends on the merger mass ratio (e.g., 
iDi Matte o et al.| 2007|). From the observationa l perspective, 
(see e.g., Ellison et al'l 1201 ll I Willett et al.l 1201 5l) . even minor 
mergers have been shown to sign ific antly contribute to lo¬ 
cal star formation dKavirail l2014aK IWillett et alJ l201.5h . In 
the case of dwarf galaxies it has been suggested that inter¬ 
actions a re responsible for t he class of blue compact dwarfs 
(BCDs) dPaudel et alJ l20L5l) : dwarf galaxies with a signif- 
icant, centrally conce ntrated young stellar population (e.g., 
iGil de Paz et alJ l2003l) . In general many dwarf systems with 
increased star formation rates have been found to be irreg- 


I 2 OI 4 I: iGarrisori-Kimimelet ^2014 ISawala et al I 2 OI 5 I) . On 

the other hand, the halo mass function is predic t ed to 
be almost completely scale -free dvan den Bosch et ^ 120051 : 
Ivan den Bosch & Jiang! 1201^ . i.e., similar for field dwarf and 
large disk galaxies. However, the subhalos of dwarf galaxies 


Ekta & Chengaluil 

2 OIOI: iLopez-SanchezI 20101; 

Holwerda et al.l 

2013; Lellietal. 

201 4bl; Knanen & CisternasI l201.5h. but in 

a number of cases no visible companion has been found 

(Brosch et alJ 20041 Ekta & Chengaluil 120101: iLonez-SanchezI 


crease in star formation and the irregular morphology are cos- 
mological gas inflows (see for example IVerbeke et al.l l2014l) 
or re-accretion of material blown out by previous starbursts, 
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and varving internal instabilities fe.s.. Meurer et al.l 

1998 

van Zee et alJ 

2001; 

Lelli et alJ l2014al: lElmegreen et alJ 

2012 

Bekki & Freeman |2002). 


In view of the above discussion, it appears plausible that 
some of these starbursts could be induce d by interactions with 
dark sa tellites. Following the brst paper bv iStarkenburg & Helmil 
( I 2 OI 5 I) here we focus on the effect of a m erger on a gas rich 
dwarf galaxy. [Starkenburg & Helmil d2015h ran a suite of collis- 
sionless simulations of dwarf galaxies and their dark satellites, 
and showed that these could severely alter the morphology and 
kinematics of the dwarf. In this paper we extend these simula¬ 
tions to dwarf galaxies with varying gas fractions and study the 
effects of the merger on the star formation rates, gas and stellar 
morphology and kinematics. The initial conditions of the dwarf 
galaxies and their satellites and the parameters for the simula¬ 
tions are described in Sect. |2] We report our results and their 
dependence on properties of the systems and the interactions in 
Sect. |3] We com pare these SPH-result s with the collissionless 
simulations from istarkenburg & Helmil (12015h in Sect. |4] and to 
observational results in Sect. |5] We conclude by discussing our 
results in Sects. |6] and give a summary of our main hndings in 
Sect.lT] 


2. Models 


We perform a suite of controlled simulations of isolated dwarf 
galaxies and mergers with their (dark) satellites. The sim¬ 
ulations are run using the OWLS version (last described 
|201(^ of the N-body/SPH-cod e Gadget- 
Springel et al.l 120011: ISpringell l2005l) with im- 
ple mentations for star form a tion a nd feedback as described 
[Schave & Dalla Vecchial (l2008h : iDalla Vecchia & Schavel 


in ISchave et~^ 
3 (based 


on 


( 120081 ) . 


2.1. Initial conditions 

The setup of the 


Springel et al.l (2005) 


in itial conditions is b ased _ on 

Schave & Dalla Vecchial (l2008l) : 
Dalla Vecchia & Schavd ( 2008) and a more com plete de¬ 
scription can be found in IStarkenburg & Helmil (l2015l) . Here we 
briefly describe the initial structure of the dwarf galaxy and of 
the satellite. The values of the structural and orbital parameters 
for all the simulations are listed in Tables [1] and 


2.1.1. Main ((disky dwarf) gaiaxies 

The host dwarf galaxies have several components including a 
dark halo, a stellar disk and a gaseous disk. T he dark matter halo 
follows a Hernquist profile (lHernauisdll990l) . 

where the parameter s po and a are set by an equivalent NFW 
(iNavarro et al.l Il996l) profile, such that the total mass of the 
Hernquist halo equals the virial mass of the NFW halo and 
their profiles have s imilar inner densities (gpo^nfw = apo.n) 
(ISpringel et alfcOOSh . 

We consider three different concentrations Cho.st = Gir/GNPw 
for the halos: 5, 9 and 15. The latter two are consistent with 
the mass - concentration relation found in large cosmolog ical 
simulations (iMaccio et al.ll2008t iMunoz-Cuartas et al.ll201 ll) for 
the mass-scale considered here. The lowest Chost corresponds to 
the best-fit NFW model for the Fornax dwarf spheroidal galaxy 


from iBreddels & Helmil (l2013l) . The velocities of the halo parti¬ 
cles are set using the distribution function of a Hernquist halo. 
Since the contribution of the disk is hereby initially neglected, 
the halo shows a slight adiabatic contraction at the start of the 
simulations, which stabilizes within approximately 0.5 Gyr. 

The stellar masses of the systems follow generally (ex¬ 
trapolated) stellar mass - halo mass relations available 
in the literature dBeh roozieLaLl 20131: Moster et al.l 120131 : 


iGarrison-Kimmel et al. 120141 : Sawala et al.H2015i) . The gas frac¬ 
tion, fg = Mgas/(M* -H Mgas), ranges from f„ - 0.3 to f„ - 0.9 
in agreement with o bservational estimates (IHuang et al.ll2012t 
iMcOuinn et al.ll20l5l) . This means that the baryonic mass can be 
quite high, and that in the most gas-rich systems Mgas = 9M*. 

Both the stellar and gaseous disk follow an exponential sur¬ 
face density profile with radius. The scale le ngths of the stel ¬ 
lar disks are close t o the values expected from iMo et al.l (Il998h : 
ISpringel & White! (Il999h for the - 1.4 x IO'^Mq and 
Chost =15 systems, and also for the lowest mass dwarf (Mvir = 
0.55 X 10^°Mo) with Chost = 5 dark matter halo, assuming that 
the disk angular momentum fration equals the disk mass frac¬ 
tion {jd - nid). For completeness we explore for two systems 
(Mvir = 0.97 X IO'^Mq and Myir = 1.4 x IO'^Mq) smaller disk 


scale lengths. For all disks the exponential radial profiles and the 
scale lengths are stable during evolution in isolation. The stellar 
disk is further described by an isothermal vertical distribution 
with a constant scale height such that O.lRd < Zo < O.SRd'. 




Md.* 

^nRlzo 


exp 



( 2 ) 


We consider thicker disks for the lowest mass simu- 

lated galaxies. This is _ in concordance with observa tions 

(lYoachim & Dalcantonll20d^ ISanchez-Janssen et al.ll2010l) and 

also expected since gas cooling is less efficient in sm aller halos 
(iKaufrnann et al.l[2007t [Robertson & Kravtso^l2008l) . 

We explore two different scale lengths of the gaseous disk in 
comparison to the stellar disk, namely Rg - Rj, and Rg - 2Rd. 
The vertical distribution of the gas is determined by requiring 
hydrostatic equilibrium. 


1 dP 50 
Pg dz dz 


(3) 


Using an effective equation of state (EOS), with yEOS.eff = 
{d In P) / (d In p), this can be written as 

^Ps pI . 4 ^ 

dz rEOS.effF dz ' 

In our simulations we use th e effective equation of s tate o f 
the multiphase ISM model bv jSchave & Dalla Vecchi^ (l2008l) : 
iDalla Vecchia & Schavel (l2008l) . The gas is governed bv the ef¬ 
fective equation of state when the density of the gas is above the 
threshold for star formation as set by the subgrid physics dis¬ 
cussed in Sect. 12.31 while at lower densities the gas is assumed 
to follow an isothermal equation of state. 

The vertical distribution of the gas can be iteratively deter¬ 
mined as a fu nction of radius usin g a fine logarithmic grid in 
the R-z plane (ISpringel et al.ll2005l) . where the surface density 
of the gas, 1.d,gas (R), is set by the exponential radial profile and 
the chosen scale length. 


( 5 ) 
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To cal culate the potential in Eq. (|4|i we follow ISpringel et alJ 
(l2005h in using a tree code and a discretized mass distribution to 
represent the disk components, adding contribution of the dark 
matter halo analytically, and use a grid of 4096 x 64 x 128 in the 
radial, vertical and azimuthal directions. 

In IStarkenburg & Helmil(l2015h we adapted the way in which 
the velocity structure of the stellar disks is set up. Contrary to 
large galaxies, for lower mass systems the epicyclic approxi¬ 
mation breaks down in a large part of the disk which, if used, 
can cause a non-physical streaming velocity. Therefore we set 
up the velocity structure only assuming the epicyclic approx¬ 
imation in the part of the disk where it is still valid and fit a 
smooth quadratical curve to the rotation curve closer in. More¬ 
over, the potential of the model dwarf galaxies is dominated by 
their dark matter halo at all radii, so we assume that the velocity 
ellipsoid is aligned with the spherical coordinate system instead 
of the cylindrical coordinate system of the disks. For the gas, the 
initial azimuthal streaming velocity is determined through the 
gravitional, pressure, and centrifugal forces; 


2 


(6) 


2.1.2. Satellites 

The satellite consists solely of a dark matter halo with a mass of 
20 % of the dwarf galaxy’s halo mass in most of our simulations. 
Its concentration is determined by following the m ass - concen¬ 
tration relation from lMunoz-Cuartas et al.l ( 201 ih . In Sect. 13.2.^ 
we explore the effect of a two or four times smaller mass and a 
higher concentration for the satellite. 

Most of the satellite orbits we consider are close to com¬ 
pletely radial with initially none or a small (prograde) tangential 
velocity, although we also explore more circular orbits (see Ta- 
ble|2]and Sect. l3.2.TT i. At the start of the simulations the satellite 
is placed at a distance of 0.67rvir from the center of the host with 
a radial velocity that is small with respect to the local circular 
velocity (so the satellite is initially close to apocenter). The or¬ 
bit is either in the plane of the host disk or has an inclination of 
30 (or in o ne case 60) degrees, similar to other studies o f disk 
thickening (IVillalobos & Helrnill2008t iMoster et al.ll2010)) . 


2.2. Numerical parameters 

We consider two different numerical setups for the host dwarf 
galaxies. For the disk3-gas (models C4 and C5) and FNX-gas 
(models El and E 2) systems that we compare t o the collission- 
less simulations of IStarkenburg & Helmil (l20L5h . we use 5x10^ 
particles to represent the dwarf’s dark matter halo and 10^ par¬ 
ticles for the gaseous and for the stellar disk irrespective of the 
value of fg. For all other hosts the 2 x 10^ baryonic particles 
are divided among gas and stars according to fg so that all gas 
and star particles have initially equal mass, and the dark halo is 
represented by 10^ particles. 

The satellite is in all cases represented by 5 x 10^ par¬ 
ticles which gives a dark matter particle mass for the satel¬ 
lite of 2 X 10^ Mq < OTsat ^ 2 X 10"^ Mq. The softening 
lengths used are 0.025 kpc and 0.016 kpc for the host halo and 
satellite respectively, and 0.008 kpc for both the gas and stars . 
The se values are chosen follo wing IVillalobos & Helmil (l2008h 
and I Athanassoula et alJ (l2000t) and produce stable systems. All 
numerical parameters are summarized in Tabled 


Table 1. Structural parameters for the host dwarf galaxies. 


Model 

Mvir 

10'°Mo 

Gir 

kpc 

c 

M* 

lO^Mo 

fid 

kpc 

O 1 


A 

5.6 

77 

9 

1.4 

0.93 

0.1 

0.5 

B 

2.2 

56 

15 

0.27 

0.78 

0.2 

0.75 

Cl 

1.4 

48 

15 

0.11 

0.78 

0.3 

0.9 

C2 





0.39 



C3 



5 


0.78 



C4 

1.4“ 


15 

0.55 


0.2 

0.5 

C5 

a 



0.77 



0.3 

D1 

0.97 

42 

5 

0.044 

0.95 

0.3 

0.9 

D2 



15 





D3 



5 



0.5 


D4 



5 


0.48 

0.3 


D5 



5 



0.5 


El 

0.55* 

27 

5 

0.22 

0.95 


0.5 

E2 

h 



0.31 



0.3 


Notes. ... denotes that the value is equal to that reported in the row 
above. 

These systems are e quivalent to the diskS-systems from 
IStarkenburg & Helmil (1201 5h but include gas. These systems are 
equiva lent to the FNX-analog-systems from IStarkenburg & Helmil 
( 1201 5h but include gas. 

Table 2. Parameters for the satellites and their orbits 


4^.sat/4^vir„ai„ “^sat 

Vr 

Vt 

’■apo/Vri“ inclination 

^^irmain 


0.2 

15 

-0.08 

0.06 

~40 

30 





~40 

60 



0 


~40 

0 


25 

-0.08 


~40 

30 



0 


~ 35 

0 




0.86 

2 





0.5 

6 


0.1 

16 

0 

0.06 

~45 

0 


25 



~45 


0.05 

17 



~40 



25 



~ 35 



Notes. ... denotes that the value is equal to that reported in the row 
above. 

The apo-to-peri ratio is defined for the first pericentric passage. For 
very radial orbits, this ratio is uncertain (by ~ 20%) because of the 
dependency on the time-sampling around this passage. 


2.3. Star formation and feedback prescription 

We use the star formation and ste llar fe edback prescription from 
ISchave & Dalla Vecchial (l2008l) and iDalla Vecchia & Schavd 
(E008l) . The star formation prescription is based on em- 
pirical laws and has very good numerical properties. 
ISchave & Dalla Vecchial (l2008h have show n how the gas 
surface densities in the Kennicutt-Schmidt law dKennicu tlll998h 
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Table 3. Numerical parameters 


Model 

-^sat 

Ndm 

xlO® 

-^gas 

xlO^ 

Ni, 

xlO^ 

moM 

lO^Mo 

^bar 

lO^Mo 

44vil-main 

A 

0.2 

1 

1 

1 

5.60 

13.8 


0.1 







0.05 






B 

0.2 


1.5 

0.5 

2.12 

5.34 

Cl -C3 



1.8 

0.2 

1.45 

5.55 

D 





0.97 

2.21 

C4 


5 

1 

1 

0.27 

5.51 

C5 






7.71“ 

El 





0.17 

3.35 

E2 






4.69'' 


Notes. ... denotes that the value is equal to that reported in the row 
above. 

(o) Pqj (Pe system C5 m^, and the value reported gives m*. = 

3.31 X 10^ Mq. For the system E2 wi* mgas and the value reported 
gives m*. = 2.01 x 10^ Mq. 


can be related to pressure by assuming that the scale height of 
the gas disk is of order of the local Jeans scale and that the gas 
is in local hydrostatic equilibrium and self-gravitating. This is 
combined with a polytropic effective equation of state for the 
multiphase interstellar medium with a slope of 7 eof = 4/3, 
which always ensures a constant Je ans mas s (so i ndependent 
of the local gas density). We follow ISchavel (l2004t) in using a 
density threshold rin - 0.1 cm“^ which corresponds to a surface 
density th reshold of ~ 10 Mp pc~^ and a temperature threshold 
of 10^ K (ISchave & Dalla VecchiallToOSh . We do include radia¬ 
tive cooling but no metals and therefore no metal-line cooling 
or chemical enrichment. 

Th e feedback model is described in iDalla Vecchia & Schavel 
(l2008l) and consists of kinetic supernova winds. The model is 
governed by the wind speed and mass loading rj - Mu,!Mi,, 
which are related to the fraction of the kinetic energy injected by 
supernovae (SNe) per solar mass, esN, also called the feedback 
efficiency as 


fu, 


Ww 

2esN 


(7) 


Following [Dalla Vecchia & Schavel (l2008h we set esN = 1.8 x 
10"^® erg which is appropriate for a Chabried (l2003h initial 
mass function, a stellar mass range 0.1 - 100 M© and all stars 
above 6 Mq ending as core-collapse SNe. The wind particles 
remain as such for t^, - 1.5 x 10^ yr and during that time are 
not able to be kicked again by another SN or participate in star 
formation, but they are not decoupled hydrodynamically. 

_ Both observational data dSchwartz & MytinI |2004 iMartir 

2005) and theoretical models ( Okamoto et al.ll2010HLagos et al. 
2013h suggest that wind velocities are lower and the mass load¬ 
ing rates are higher for low mass galaxies and low star formation 
rates. Therefore, rather than making a single choice for the nu¬ 
merical parameters, we experimented with a selection of differ¬ 
ent parameters values as described in the next section. 


2.4. Evolution in isolation 

2.4.1. Exploration of subgrid parameters 

To identify a reasonable default model for the dwarf galaxies 
we explored a small set of combinations of the large parameter 
space and varied two of the stellar feedback parameters, namely 
the wind speed (30, 90, 120, 240, and 600 km/s) and the mass 
loading T] (2, 8 , 10, 20, 30, 40, 50). 

Our goal is for a galaxy in isolation to have a fairly continu¬ 
ous star formation rate on a timescale comparable to that of the 
merger. This is important to be able to pin down to the effect of 
the merger. However since we do not model fresh gaseous infall 
all star formation rates will decline due to gas depletion. We also 
pay attention to ensure that the gas is not all converted into stars 
or blown away by stellar feedback too quickly. Nonetheless, a 
tight correlation exists between the initial gas mass and its distri¬ 
bution and the early star formation rate. Via the threshold for star 
formation, this also depends on the gravitational potential (and 
hence on the mass and concentration of the dark matter halo and 
the mass, scale length, and scale height of the stellar disk). 

The star formation rate over a period of 2 Gyr for the dwarf 
galaxy model Cl (My^ = 1.4 x lO'^’Me, /gas = 0.9, Rg - Rj, 
and zo = 0.3Rd), for varying feedback parameters is shown in 
Fig.IU The star formation rates for all prescriptions are within 
a factor of a few of each other and in all cases star formation is 
sustained over a long timescale. For the more extended gas disk 
(Rg - 2Rd, dashed dark green) the star formation rates are about 
a quarter of the Rg = R^-case (black). For the most efficient feed¬ 
back (vw - 240km/s and 77 = 30 so Fwind/^SN = 96%, dashed 
purple) the gas disk is blown apart right after the onset of star 
formation and subsequent star formation occurs at a similar rate 
as for the more extended gas disk. The periodicity seen in the 
first ~ 1 Gyr is driven by the fact that star formation starts in the 
whole disk at the same time. This is enhanced if the disk is not 
perfectly centered in the potential but is reduced at later times 
and for lower star formation rates. The amplitude of the oscil¬ 
lations in star formation also depends on the star formation and 
feedback parameters. For example the amplitude of the period¬ 
icity is larger for higher mass loading of the wind (compare for 



time [Gyr] 


Fig. 1. Star formation rates for model Cl (Myi, = 1.4 x IO'^Mq) with 
different values of the feedback parameters: wind velocity and mass 
loading depicted in different colors as indicated by the inset for Rg = Rj, 
except for dark green where Rg = 2Rj. 


Article number, page 4 ofll9l 








































T. K. Starkenburg et al.: Gas and star formation in minor mergers of dwarf galaxies with dark satellites 



Fig. 2. Face-on (top) and edge-on (bottom) density contours of the gas in the disk of model A with Rg = IRj, evolved in isolation. The contour 
levels are at 0.25, 1, 4, and 16 x 10^° N cm“^. The local star formation regions (and their amplitude) are shown in green at different moments 
during the simulation. 


example the gray, black, and blue star formation rates in Fig. [Til. 
Reassuringly, the evolution is very similar for settings with com¬ 
parable feedback efficiencies (see for example the red, yellow, 
and black curves in Fig.[T]all with Fwind/FsN ~ 24%). 

We tested the influence of decoupling the wind particles from 
star formation by comparing with a run without decoupling, and 
found that this does not change the star formation rate and evo¬ 
lution of the gas disk. We further explored the effects of using 
higher values for the density thresh old (1 and 10 cm~^ ), as 
proposed in the recent literature (e.g.. lGovernato et alJ2010l) . but 
found that this restricts star formation to the very center of the 
galaxy and does not enhance the fragmentation of the disk gas 
as we do not include metal-line cooling, and therefore does not 
lead to a more realistic system. 

As default model we take the star formation and feedback 
parameters corresponding to the black curve in Fig. [T] for Rg - 
Rd (dark green for Rg - 2Rd). The default wind velocity is = 
120 km/s with a mass loading of ?; = 30, resulting in fg, - 24%. 

2.4.2. The default model 

Figure [^ shows snapshots of the gas distribution and state for 
one of the most massive dwarf galaxies, model A (Myir = 
5.6 X 10'°Mo and M* = 1.4 x IO^Mq), evolved in isolation. 
For all plots in this paper the plane of the disk is defined as the 
plane perpendicular to the angular momentum vector of the in¬ 
ner 50% of the initial stellar disk particles. As mentioned above 
the stellar feedback parameters are set such that the star forma¬ 
tion rate and the disk of the dwarf galaxy are reasonably stable 
for several Gyr. Since the gas in the outskirts of the disk is not 
dense enough, most of the star formation takes place within the 
central ~ 2kpc(~ Rg - 2Rd for this model). Due to the stellar 



0.0 0.5 1.0 1.5 2.0 

time [Gyr] 


Fig. 3. Star formation rates (top panel) and central gas masses (bottom 
panel) for all systems in isolation with Rg = Rd for models A (orange), 
B (red). Cl (blue), D1 (black), and with Rg = 2Rd for models A (purple) 
and Cl (green). Note that although model Cl (blue) is less massive than 
B (red), the star formation rate is higher for Cl because it has a higher 
gas mass in the center due to its higher initial gas fraction. 

feedback which blows gas away that eventually falls back again, 
the star formation is patchy and locally bursty. 
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Fig. 4. Evolution of model A (M* = 1.4 x 10* Mq) with Rg = 2R^ merging with a 20% mass satellite on a co-planar, very radial orbit. The top 
rows show the face-on and edge-on, respectively, view of the gas in the disk (contours at 0.25, 1, 4, and 16 x 10^® N cm“^) with the gas that is 
currently forming stars highlighted in green (see colorbar for relative values). The bottom panels show the old stellar component in red, and newly 
formed star particles in black, at different times during the merger. 





gas edgeon 




stars edgeon 


The star formation rates (SFR) for a subset of the dwarf 
galaxies run in isolation using our default star formation and 
feedback parameters are plotted in the top panel of Fig. |3] The 
average amplitude of the star formation rate clearly depends on 
the total mass of the galaxy, the amount of gas in the disk and its 
initial extent. Indeed, for a more extended disk set up in equilib¬ 
rium, a smaller fraction of gas will be above the density threshold 
for star formation. 

The amount of gas in the central part of the disk is shown 
in the bottom panel of Fig. [3] That both the gas fraction and 
the mass of the dwarf galaxy are important is made explicit by 
comparing model B (Myir = 2.2 x M* = 2.7 x IO’Mq 

and fg = 0.75, red) and model Cl (Myir = 1.4 x lO'^’Me, M* = 
1.1 X 10^ Mq, fg - 0.9, blue). The bottom panel of Fig. |3] shows 


that although the more massive dwarf galaxy model B has 2.5 
times more mass in dark matter and stars initially, the lower mass 
system model Cl actually has a slighty higher gas mass (Mgas = 
9.9 X 10^Mq versus Mgas = 8.2 x \QpM q, due to the higher gas 
fraction) and a higher star formation rate. After 2 Gyr the stellar 
masses of the two systems are similar and also the gas fractions 
are more comparable than initially. 

Although the simulated systems have slowly declining star 
formation rates in general (as a result of the lack of exter¬ 
nal gas infall), they ag ree quite well with observa t ions, also in 
terms of gas content (iHuriter^&Elmegrera J200^Weisze^^ 
20121 lOtt et all IMl l^rren et al.1 1201 2t [Huang et al.l 120121: 


McOuinn et al.ll2015l) . However, the lower mass disks tend to be 


more extended, slightly thinner, and have lower surface bright- 
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Same as Fig.|4]for model D1 (M* = 4.4 x 10 ^Mq with Rg = Rj), merging with a 20% mass satellite on a 30-degrees inclined, very radial 


ness compared to observations. We make a more detailed com¬ 
parison with observations of the isolated and merger remnant 
dwarf galaxies in Sect. |5] 


3. Starbursting dwarfs as the result of an interaction 

In this Section we compare first the behaviour of one of the high¬ 
est (model A) and one of the lowest (model Dl) mass dwarfs, as 
they interact with a dark satellite. We describe the general prop¬ 
erties and then address in more detail how the interacting sys¬ 
tems vary depending on the properties of the satellite, its orbit, 
and the host. In all cases explored the evolution is quite differ¬ 
ent from the dwarfs in isolation, as the majority of the systems 
experience starbursts of varying strength. 


3.1. Two examples 

Figure |4] shows a series of snapshots of the evolution of the sys¬ 
tem of Fig. model A (M* = 1.4 x IO^Mq) with Rg - IRd, 
now merging with a satellite on a co-planar very radial orbit. 
The stellar disk thickens, tilts, and develops tidal arms and rings. 
The gaseous disk depicts minor tidal arms and becomes quite 
asymmetric. Most importantly gas is driven into the center of 
the dwarf galaxy due to tidal torques, which leads to a strong 
increase of the star formation rate as we discuss below. 

Fig. |5] shows the evolution for the lower mass model Dl 
(M* = 4.4 X IO^Mq) merging with the satellite on a very ra- 
dial 30 degrees inclined orb it. Similar to what was found in 
IStarkenburg & Helmil (1201 5h . for this smaller mass object the 
effect of the satellite on the stellar disk is much stronger. The 
stellar morphology of the remnant is spheroidal, for both the old 


Article number, page 7 ofll9l 


























A&A proofs: manuscript no. MergersDwarfsandDarkSatellitesGasandSF 




Fig. 7. Face-on (top) and (bottom) edge-on gas contours for mode 
A with Rg = 2R^ merging with a 20% mass satellite on a radial, co- 
planar orbit depicted in Fig. |4l initially (far left) and after 1 (middlf 
left), 2 Gyr (middel right) and 3 Gyr (far right). The plane of the disl 
in determined by the angular momentum of the inner 50% of the initia 
stellar particles, so the rotation in the gas can be in a different plane. 


Fig. 6. Star formation rates with respect to the average of the SFR in 
isolation for the first Gyr, for the systems A1 (see Fig. 0 and D1 (see 
Fig. [Hi experiencing a 20% mass ratio merger. A number of pericentric 
passages of the satellite during the different simulations are indicated 
by the arrows at the top of the figure. 


and the newly formed stellar populations (see the rightmost pan¬ 
els of Fig. |5]l. The gas disk is severely disturbed, although no 
strong tidal tails form. As can be seen in the third panel of the 
top two rows, the presence of the satellite can cause off-center 
starbursts. 

Fig. |6] shows clearly the effect of the merger on the SFR for 
these two systems. Both for model A (blue) and D1 (black) the 
star formation rates are increased by factors ~ 3 up to 12. This 
happens not only during the merger process itself with sharp 
strong peaks at or just after pericenter passages, but more strink- 
ingly for a rather extended period. The relative amplitude of the 
increase in the SFR is similar for both systems, although in an 
absolute sense, the more massive object naturally has a higher 
SFR, reaching values of 0.04M© yr“*. Note as well from Figs.|4] 
and |5] that in both cases the gas disk at the end of the simulation 
is much more compact than initially. 

In Figures|2]andOwe plot the evolution in the rotational ve¬ 
locity for the two systems shown in Figs. |4] and |5] respectively. 
The rotational motion of the gas disk of model A seems to in¬ 
crease first (see the left middle panel of Fig.|7]i and then strongly 
decreases as the gas disk shrinks until the significantly smaller 
maximum circular velocity ~ 30 km/s is reached. Nonethe¬ 
less, the velocity held generally remains rather robust during the 
merger. For the lower mass system, which depicts a smaller am¬ 
plitude of rotation initially, the velocity held is much less con¬ 
spicuous and ordered during and after the merger. Note that at 
the end, the maximum rotation signal is found for what we have 
dehned as the “face-on” view of the system, although this char¬ 
acterization is debatable given the spheroidal shape of the rem¬ 
nant. 

3.2. Variation in properties of the sateiiite 

For a signihcant effect, the satellite must reach the stellar and gas 
disks of the dwarf galaxy within a short timescale to be apparent 
in our simulations. Satellites on orbits that are close to circular 
will take a longer time to sink to the center. This is longer than 
the few Gyr run-time of our non-cosmological simulations to 
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Fig. 8. Same as Fig.|7]now for the D1 system depicted merging with a 
satellite as shown in Fig. [5] 


limit environmental and cosmological effects, such as the lack of 
cosmic gas inflow. Moreover, the strength of the perturbation de¬ 
pends on the average density ratios of the satellite to the host. To 
explore the dependencies of the mergers on the satellite’s proper¬ 
ties we consider a number of different orbits, and satellites with 
different concentrations (and central densities), for three differ¬ 
ent mass ratios: 1 : 5 (our default), 1 : 10, and 1 : 20. 

3.2.1. The satellite orbit 

Figure |9] compares the star formation rates for the dwarf galaxy 
model A experiencing a 1:5 merger with the same (high concen¬ 
tration) satellite on three different orbits. For less radial orbits the 
disk forms large tidal spiral arms (compare for example Figs. IA.il 
and lA.2l in the Appendix) and the subsequent increase in star for¬ 
mation only starts when the satellites comes within ~ 3 kpc of 
the center. For the most extreme example shown here (the red 
line in Fig. |3 this happens only at 3.5 - 4 Gyr, i.e. at the end of 
the simulation run. 

For more circular orbits the spiral arms that are generated 
in the gas (and stellar) disk are more pronounced. Nevertheless, 
for the radial orbit the gas is much more concentrated after the 
merger, which suggests that more gas is funnelled to the center 
due to tidal torques. 

We also experimented with different inclinations for the 
satellite orbit with respect to the plane of the disk (see Figs. lA.3] 
and I A. 41 for examples). In general we find that an inclination of 
30 degrees has a larger effect on the morphology and kinemat¬ 
ics of the stellar disk (see also our collissionless simulations de- 
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Fig. 9. Star formation rates for model A with Rg = 2Rj: in isolation 
(black) and experiencing a 1:5 merger with a high concentration, c = 25, 
satellite on a planar orbit: radial with first pericenter within the stellar 
disk (blue), less radial with first pericenter just outside the stellar disk 
(light green), and even less radial with fperi ~ 20 kpc (red). A number 
pericentric passages of the satellite during the different simulations are 
indicated by the arrows at the top of the figure. 



Fig. 10. Star formation rates for model A with Rg = Rj: in isolation 
(black), and with a satellite on a planar radial orbit with = 

0.2 (Csa, = 15 blue; c^a, = 25 green), = 0.1 (c^a, = 16 red; 

Csat = 25 orange), and = 0.05 (Csat = 17 purple; c^at = 25 

gray). A number pericentric passages of the satellite during the different 
simulations are indicated by the arrows at the top of the figure. 

scribed in IStarkenburg & Helmill2015l) but drives less gas to the 
center and produces a smaller increase in star formation rates 
compared to a co-planar orbit. 

3.2.2. Satellite mass and concentration 

For our default setup the satellite mass is 20% of the virial mass 
of the host dwarf galaxy. This choice is motivated by the fact that 
interacti ons with such objec ts can be devastating and may not be 
so rare. iHelmi et"^ (1201 2h have estimated that dwarf galaxies 
experience on average 1.5 such encounters in a Hubble time, but 
encounters with smaller mass objects are certainly more com¬ 
mon in the ACDM cosmological model. 



time [Gyr] time [Gyr] 

Fig. 11. Star formation rates, relative gas fractions, relative stellar 
masses, and gas masses in the central parts for model A with Rg = Rj 
and Rg = IRj in isolation (black and blue, respectively) and during the 
merger (green and red, respectively). The arrows at the top of the top- 
left figure indicate a few pericentric passages of the satellite during the 
simulations shown. 



Fig. 12. Star formation rates for models C and D with Rg = R^, /gas = 
0.9, zq = 0.3Rd, and a dark matter halo concentration of either c = 5 (C3 
and D2) or c = 15 (Cl and Dl), during the 1:5 merger with a satellite 
on a prograde orbit with an inclination with the plane of the disk of 
0 (for models C) or 30 degrees (for models D). A number pericentric 
passages of the satellite during the different simulations are indicated 
by the arrows at the top of the figure. 


We explore the effect of 1:10 and 1:20 mergers for the most 
massive dwarf galaxy (model A) with Rg = 2Rd, and place the 
satellites on a radial, planar orbit. We cons ider satellites follow¬ 
ing th e mass-concentration relation from iMunoz-Cuartas et"^ 
(1201 Ih and also having a higher concentration of Csat = 25. 

The resulting SFR are shown in Fig. [TO] In all cases the am¬ 
plitude of the starburst depends most strongly on the mass and 
secondly on the concentration of the satellite, although the latter 
dependance is weaker the lower the mass of the satellite. Note 
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Fig. 13. Vertical (top panels) and radial (bottom panels) density profiles (solid lines) during the merger simulations with intervals of 1 Gyr for the 
FW-like system: fully collissionless (left), with /gas = 0.3 and /gas = 0.5 (middle), and with /gas = 0.5 including the newly formed stars (right). 
The dotted lines show the maximum likelyhood fitted profiles: isothermal (sech^) and exponential in the top and bottom panels, respectively. 


also that the onset of the major starburst is later for smaller mass 
satellites. This is because the satellite sinks in more slowly, i.e. 
the pericenter passages where star formation is triggered occur 
later. 

For all satellites, higher concentration leads to slower mass 
loss, hence to more damage to the host. This results in a larger 
increase in its SFR and to a stronger morphological disturbance 
of stellar and gaseous disks, with the final gas distribution being 
more centrally concentrated. While gas gets blown out of the 
disk in all cases, in the interaction with the highest mass and 
Csat = 25 satellite the gas reaches ~ Tryu, about twice as far as 
for the other cases. 

In summary, in all cases (although very minor for 
ATsat/Myir^^^ = 0.05), there is a first strong peak of star for¬ 
mation at the first pericentric passage of the satellite, and a 
more extended in time starburst, also driven by the merger. Sec¬ 
ondary peaks associated to subsequent pericentric passages are 
also present but are generally less conspicuous. 

3.3. Influence of properties of the host dwarf galaxy 

The structure of the dwarf galaxy itself can significantly alter 
the effects of the minor merger on the gas and star formation. 
We vary the halo mass and concentration, the gas fraction and 
extend of the gas disk, and the scale length and scale height of 


the stellar disk. Of these, varying stellar disk parameters causes 
only very minor differences on the gas and star formation. 


3.3.1. Gas distribution 

We explore now the evolution of a disky dwarf galaxy when the 
gas disk has different extent than the stellar disk, since this char¬ 
acteristic is often seen in large spiral galaxies. A comparison be¬ 
tween the merger properties for model A with Rg = 2Rd and with 
Rg = Rd can be seen in Fig.ITTI 

Due to the lower initial star formation rate for the more ex¬ 
tended (and hence lower surface density) gas disk, this system 
has a larger gas reservoir at the time of the merger (see top right 
panel). The bottom right panel of Fig.[TT]shows that the increase 
in gas mass near the centre is much higher for the initially more 
extended disk during the merger (red curve) than for the less 
extended gas disk (green). However, the final increase in stellar 
mass is still higher for the system with Rg - Rd, mostly due to 
its initially larger star formation rate. 

In general more extended disks have lower star formation 
rates both in isolation and during the merger. This is because the 
same amount of gas is distributed over a larger area and therefore 
the amount above the star formation threshold is much lower. 
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3.3.2. Host’s dark matter halo mass and concentration 

Figure [12] shows the star formation rates for the dwarf galax¬ 
ies models Cl and C3, and D1 and D2 which have different 
dark matter halo concentration and mass. Even though the dif¬ 
ference in halo mass between models C and D is only a factor 
1.4, due to the steepness of the halo mass - stell a r mass abun- 


dance matching relations (iBehroozi et al. 2013 

iMoster et aLl 

l2013l iGarrison-Kimmel et al. 20141 ISawala et al. 

20151) the dif- 


ference in stellar mass and gas mass is a factor 2.5. For these 
runs the stellar disk is thick initially (zo = 0.3Rd), Rg - Rd, 
/gas = 0.9, and the satellite has a mass of 20% of the host halo 
on radial orbit with an inclination of 0 or 30 degrees. 

As noted before, higher (central) mass implies initially 
higher SFR, and thus also during and after the merger (compare 
black and green, and blue and red curves). A less concentrated 
host halo initially has a puffier gas disk and lower SFR, but dur¬ 
ing a merger the enhancement in the SFR is larger, as depicted 
by the black curve in Fig. [12] for model C3. Therefore we see 
that the susceptibility of the system to a merger depends both on 
virial mass of the host and on its concentration. 



Fig. 14. Star formation rates versus stellar masses and B- magnitude in 
interva l s of 1 Gyr in compar ison t o the dwarf galax i es fromlWeisz et ^ 
(l2012h . iHuang et all (1201 2li and iMcOuinn et all (|2 01^ for the left 
panels, and i n add ition for the right panel f r om Ott et al.l ( 1201 2h . 

I Warren et all (l2012li . and iHunter & Elrnegre^ (l2006li . Note that the 
stellar mass grows during the simulations. 


4. Comparison with coilissioniess runs 

In this section we compare the dwarf galaxies models C4 and C5 
(disks-gas) and El and E2 (FNX-gas) wi th the results from the 
“equiv alent” coilissioniess simulations in IStarkenburg & Helmil 
(l20L5h . Our focus is on how the presence of the gas changes the 
characteristics of the merger remnants. 

The morphological and kinematical changes to the stellar 
component of both dwarf systems are very similar to the fully 
coilissioniess runs. We exemplify this in Fig. IA.5I of the Ap¬ 
pendix, and more quantitatively in Figure [T3] This hgure shows 
the vertical and radial density prohles of the stellar disks of the 
initiall y thick FNX-analog system from IStarkenburg & Helmil 
(1201 5h (left panels) and for the counterparts with fg - 0.3 and 
fg = 0.5 (central panels). In all three cases the strongest changes 
occur during the first pericenter passages of the satellite, within 
2 Gyr of the start of the simulations. The differences in the struc¬ 
tural evolution of the stellar disks are surprisingly small. The 
major difference is that for the fully collisionless system the ra¬ 
dial profile of the disk is slightly more unstable and has more 
substructure. Even the newly formed stars in the runs including 
gas do not significantly change the stellar density profiles as can 
be seen in the rightmost panels of Fig. but they give rise to 
a small central bulge (compare the bottom panels). For the more 
massive diskS (models C4 and C5) the inclusion of gas also has 
a small effect, and leads to a slightly thicker stellar disk in the 
post-merger phase, more so for higher gas fractions. 

The effect of the merger on the kinematics of the stars do not 
differ much when gas is included in the disk. For example, for 
the FAW-like system in both cases the rotation decreases over 
the whole disk while the circular velocity increases due to the 
accretion of the dark matter satellite. All velocity dispersions in¬ 
crease, and even increase slighly more in the presence of a gas 
disk. This might be due to the fact that the gas disk is not very 
thin due to the low mass of the dwarf galaxy halos and there¬ 
fore exerts less of a pull toward the disk midplane compared to 
a purely stellar disk. 

Therefore, the presence of gas in the disk has no significant 
influence on the effects of minor mergers on dwarf galaxies. This 
is in striking contrast to what happens in similar simulations for 
higher (Milky Way) mass galaxies (see discussion in Sect.[6]l. 



Fig. 15. Gas fractions versus stellar mass es in interva l s of 1 Gyr 
in comparison to the dwarf galaxies from iHuang et all (l2012h . and 
iMcOuinn et all ( 120151) . Since the stellar mass grows during the simu¬ 
lations and that we do not include fresh gas inflow, the gas fraction of 
our simulated systems necessarily decreases. 

5. Comparison to observations 

We now compare the properties of our simulated dwarf galax¬ 
ies, in isolation and during the merger, with the observational 
samples described in Table [4] This comparison includes systems 
from Hl-selected, as well as mass-selected samples and general 
compilations. 

The left panel of Fig. [14] shows the SFR versus stellar mass 
Mt for the systems disc ussed in this paper c ompared to observa ¬ 
tion al values reported infWeisz et al.|il2012h . lHuang et ^(l2012h . 
and IMcOuinn et all (1201^ . For each simulation the evolution of 
the SFR and M, are plotted at four epochs in steps of 1 Gyr, 
from 1 Gyr after the start of the simulation. The black symbols 
correspond to the runs in isolation while in color we show the 
mergers. Overall the star formation rates match very well those 
from observations for their stellar masses. 

To estimate the luminosities and surface brightness of our 
simulated galaxies we need to assume mass-to-light ratios for the 
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Fig. 16. Central surface brightness (top) and disk scale lengths (bot¬ 
tom) versus absolute magnitudes in intervals of 1 Gyr compared to ir¬ 
regular dwarf galaxies (black), blue compact dwarfs (dark blue), and 
Magellanic spirals (gray) from iHunter & Elmeereenl l l2006h . Note that 
the stellar mass grows (mostly in the center) and so the magnitude and 
disk scale length decreases and the central surface brightness increases 
during the simulations. 


stellar particles. One possibility for the stellar particles formed 
during the simulation would be to use stellar population models 
as we know their ages. For the stellar particles initially present 
in the simulation however the age distribution would be arbi¬ 
trary. Therefore, we instead fitted a linear relation to the log M*- 
B magnitude distribu tion for the galaxies in the observations 
of IWeisz et^ (1201 2h and iMcOuinn et alJ d20L5h and apply it 
to the simulated dwarf galaxies. The resulting estimated val¬ 
ues are roughly consistent with an average mass-to-light ratio 
of Mi,!L b,* ~ 0.5, which is quite reasonable as the light is very 
strongly dominated by young stars. The right panel of Fig. [14] 
shows SFR against the B-magnitude for our systems and for a 
larger sample of observations. Good agreement is also found in 
this case. 

Figure [15] shows that also the initial gas fractions and their 
evolution in time agrees very well with observations. The down¬ 
ward trend seen in the simulated systems is due to the fact that 
we do not include fresh gas inflow. However, sampled at an ar¬ 
bitrary point in time, the match is quite remarkable. 

We compute radial surface brightness profiles for the face- 
on disks and fit an exponential profile to derive the disk scale 
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Fig. 17. Distribution of “projected” axis ratios and rotational velocity 
for 100 random inclinations at intervals separated by 1 Gyr during each 
of the simulation runs, compared to galaxies with 
100 k m/s from the catalog of nearby galaxies dKarachentsev et afl 
120041) (gray). Darker colors indicare a higher density of objects in this 
plane. 


length and half-light radius . This process is similar to that of 
iHunter & Elme green] (l2006b and we compare our disks with 
their observational results. We derive central surface brightness 
using the mass/luminosity within the innermost, 0.1 kpc, bin in 
projection. The results are shown in Fig. [16] It is clear that while 
the higher mass disks match very well with the observations, the 
lower mass systems seem to be more extended and fainter, com¬ 
pared to the observations. There are two main possible reasons 
for this: low surface brightness, extended systems are harder to 
observe. On the other hand, the method we used to compute the 
initial size of the systems and to set up our simul ations (using the 
disk mass and properties of the halo following iMo et aHll998h 
might break down for lower mass systems, as they typically are 
thicker and deviate from being thin disks. To estimate the effect 
of this we have also run the two lower mass systems with initial 
disk scale lengths half their default values and, not surprisingly, 
we find that they agree better. 

In Fig. [17] we compare our “observed” axis ratios and rota¬ 
tional velocities to the observed axis ratios of all galaxies with 
wsQ < 200 km/s from the catal og of nearby galaxie s (450 when 
the Milky Way is excluded) bv iKarachentsev et al.l (l2004l) . The 
intrinsic axis ratios are computed by fitting an ellipsoid to the 
stellar particle distribution (normalized b y the ellipsoidal dis - 
tance of the particles within this ellipsoid. lAllgood et al.ll20()^ . 
The ellipsoids are centered around the center of mass defined 
by the stellar component and the axis ratios computed for parti¬ 
cles within major axis length of a ~ 1.4 kpc, but we find similar 
results when considering particles within a ~ A kpc and even 
a ~ 14 kpc. These intrinsic minor-to-major axis ratios are in 
agreement for most of the merger remnants, but the default ini¬ 
tial disks are often thinner than the observation^ However, most 
systems will not be observed edge-on, and therefore to provide 
more realistic estimates of “observed” axis ratios we put each 


* Initially thicker disks (zo = 0.5Rrf), which might be expected for 
dwarf galaxi es that form in gas that has cooled less efficient tha n in 
larger disks dKaufmann et al.Tl2007l : [Robertson & Kravtsov! [20081) . al¬ 
ready agree well with the observational points. 
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Table 4. Observational samples used 


Reference 

Sample 

Properties used 

Hunter & Elmegreen (20041 and 
Hunter & Elmegreen (20061 

94 Im, 24 BCD, and 18 Sm 

Mm, My, B -V, Rd, SER, po 

Weisz et al. (20121 

185 galaxies from 

the Soitzer EVE survev (e.g.. Dale et al. 200911 
and 1IHUGS (Kennicutt et al. 2008: Lee et al. 20111 

M*, Mb, SER 

Huang et al. (20121 

229 low HI mass galaxies from 

the ALEALEA survev (Giovanelli et al. 20051 

M*, Mm, SER 

Ott et al. (20121 

VLA-ANGST survey: 35 galaxies from 
the ANGST survev (Dalcanton et al. 20091 

Mb, SER 

Warren_e^^ (20121 

31 nearby low-mass galaxies taken from 

THINGS (Walter et al. 20081 and 

the VLA-ANGST survev (Ott et al. 20121 

Mb, SER 

McOuinn et al. (20151 

12 galaxies from the SHIELD survev (Cannon et al. 20111 

Mm, M*, SER, Mb 

Karachentsev et al. (20041 

all-sky catalog of basic optical and HI properties of 
(451) neighboring galaxies with D < 10 Mpc 
or radial velocity Vlg < 550 km/s 

Vrot, bja 


system at 100 random inclinations and compute projected axis 
ratios at five evenly spaced epochs during the simulation. 

For the rotational velocities, we assume that for the observa¬ 
tional sample Uiot ~ 0.5 x w^o and compute the rotational veloci¬ 
ties for our modelled systems as the mean between the minimum 
and maximum value of the velocity maps (as in Fig. |7]l, which 
give the velocities in the disk within the column density contour 
of 2.5 X lO'® N cm“^ (which is close typically to ~ 2Rd). 

Inspection of Figure [T7] shows when including the (random) 
inclination angle, the sample of simulated dwarf galaxies agrees 
very well with the observed sample. Surprisingly, even the trend 
that the lower-mass stellar disks are more extremely perturbed 
by the minor merger is apparent. 

Overall, our systems agree well with properties of observed 
dwarf galaxies. 

6. Discussion 

In this paper we show the effects of minor mergers with a dark 
satellite on gas-rich dwarf galaxies. The isolated systems a care¬ 
fully set up and we have explored the effects of varying the sub¬ 
grid parameters on the results. Our star formation and feedback 
prescriptions lead to reasonable dwarf galaxies when evolved in 
isolation, and their properties are fairly stable against varying 
these parameters. Only extremely efficient feedback gives a com¬ 
pletely different and unrealistic evolution. We have also carried 
some of the merger experiments for different feedback schemes 
and found that the effects of the merger on the gas and star for¬ 
mation evolution to be very robust against such changes. 

When comparing the fully hydrodyna mic models to the col- 
lission less minor merger simulations of IStarkenburg & Helmil 
(1201 5h we have found that the evolution of the stellar compo¬ 
nents of the dwarf galaxies is rather similar. Even though the 
gas absorbs and dissipates some of the energy injected via the 
merger, the morphological transformations that the stellar disks 
experience are still very important. This is in contrast to what 
has been reported in the literature for Milky Way-size systems 
wher e the effect is significantly reduced whe n gas is included 
te.g.- lHonkins et ai]l2009t iMosto et al.ll201fi) . This is probably 


due to the gas being much colder in larger disk galaxies and 
therefore having a stronger stabilizing effect on the stars in the 
disk. 

An important question is how often the process discussed 
in this paper would happen for dwarf galaxies in different en¬ 
vironments, at different redshifts and with different masses. As 
the CDM halo mass function is almost s cale-free (with small 
differences due to halo formation times Ivan den Bosch et ^ 
120051 : Ivan den Bosch & Jiangl l2014ll . dwarf galaxies will have 
a spectrum of pert urbers very similar to that of an L* galaxy. 
iHelmi et"^ ( 1201 2|) estimate that the number of minor mergers 
for low galaxy efficiency systems (M£;/(Mvii- x/bai) = 5%) where 
the satellite is at least as massive as the disk at pericenter, with 
the pericenter within 30% of the virial radius, is ~ 1.5 within a 
Hubble time. The simulations discussed in this paper were for 
systems with even smaller galaxy efficiencies experiencing 1:5 
and 1:10 mergers with dark satellites, albeit on very radial or¬ 
bits. This means that almost every dwarf galaxy should have 
experienced a minor merger with major effects during its life¬ 
time. Predictions and their dependencies on environment, red- 
shift and dwarf galaxy mass will be presented in a forthcoming 
paper (Starkenburg et al. in prep.). 

Several defini t ions o f a starburst exis t in the literature (e.g., 
Knanen & .lamesi l2009t iMcOuinn et al.l l2010l : iBergvall et al.l 
2015 1. The birthrate parameter, f>=SFR/(SFR) is often used, and 
compares the current SFR or the peak of the burst to the average 
SFR over the past Gyr or even the lifetime of the system. For the 
more massive systems in our simulations the SFR in isolation 
is very similar to their initial stellar mass divided by a Hubble 
time, and the birthrate goes up to ~ 3, or even b ~ 10, dur¬ 
ing the merger, depending on the configuration. The lower mass 
systems have typically more bursty star formation histories even 
in isolation. We may characterize b ~ \ between these small 
bursts, b ~a few during bursts in isolation, but during the merger 
we find > 10. Therefore the increase in star formation rates 
the simulated systems experience can be qualified as starburst 
events. Also when computing the gas consumption timescale, 
Tgas = Mgas/SFR, there is agreement with the literature. For the 
star formation rates in isolation, the gas consumption timescales 
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are long >10 Gyr, while these drop significantly by fac¬ 

tors of a few during the merger events. 

Although observations suggest that only a small fraction of 
dwarf gal axies are curren tly experiencing a starburst (~ 6% ac¬ 
cording to lLee et ^l2009h . this might be just a fraction of the to¬ 
tal numb er of dwarf galax ies that experienced starbursts in their 
lifetimes (iLee et al.ll2009l) . 


7. Conclusions 

We performed a suite of controlled, minor merger simulations 
between carefully set-up gaseous dwarf galaxies and their (dark) 
satellites. These interactions can give rise to a strong increase 
in the star formation rates in the dwarf galaxies. The increase 
is in the form of large sharp bursts during pericenter passages, 
as well as of extended boosts due to tidal torques funnelling gas 
toward the center. The gas and stellar disks show severely dis¬ 
turbed morphologies in most cases, especially for lower mass 
hosts experiencing a 1:5 merger. The gas disks can develop grand 
tidal tails and their remnants depict a much more concentrated 
final distribution in some cases. For the lowest mass systems ex¬ 
plored (M, < 1.1 X 10^ Mq), the merger can completely destroy 
the stellar disk. These objects become spheroidal-like and have 
bursty star formation in their center. In contrast to simulations 
of Milky Way-like systems, the presence of gas in the disk of 
the dwarf galaxies does not diminish the effect of the merger 
on the stellar component. Our simulations that include gas show 
that the strong heating and evolution of the stellar disk is almost 
completely identical to the collisionless case. 

We have explored the dependence of the mergers and their 
remnants on the host and satellite masses, dark matter halo con¬ 
centrations, satellite orbits, gas fractions, and structure of the 
stellar and gas disks. To have a significant impact, a dark satellite 
must have at least 10% of the mass of the host. More generally, 
the strength of the merger’s effects decreases with lower satel¬ 
lite masses, but this depends on satellite halo concentration; a 
very dense low-mass dark satellite can survive longer than one of 
lower concentration, and therefore have stronger effects. Orbits 
in the plane of the disk cause stronger starbursts, while inclined 
orbits perturb the stellar components more effectively. Satellites 
on radial orbits cause stronger starbursts than those on more cir¬ 
cular ones. 

When the hosts have lower concentrations, the 
merger induces stronger morphological changes (as in 
IStarkenburg & Helmil l2Q15li but in general also lower star 
formation rates. This is because SFR correlates with gas density, 
which depends in turn on the total mass distribution in the region 
probed by the gas disk. Also the gas fraction and distribution, 
as well as stellar disk masses and distributions, determine the 
amount of gas that has high enough densities for star formation, 
thereby directly affecting the amplitude of the starburst that a 
dwarf experiences during a merger. 

Both our initial systems and their remnants compare well 
with the observational properties of a large selection of irregular 
dwarf galaxies and blue compact dwarfs. Even systems that are 
strongly perturbed as a result of a merger with a dark satellite fall 
within the scatter seen in the observations. This implies that such 
events might well be happening but may not be fully evident. We 
have yet to identify the “smoking gun” of the dark merger sce¬ 
nario. However, this also shows that the effects of interactions 
with dark satellites, which are naturally expected within a CDM 
cosmology, are likely to play a role in the diversity of the dwarf 
galaxy population. 
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Fig. A.l. Evolution of model A with Rg = 2Rj, merging with the c = 25 satellite on a planar, very radial orbit, at different times during the merger. 
The top rows show the face-on and edge-on, respectively, view of the gas in the disk (contours at 0.25, 1, 4, and 16 x 10^® N cm“^) with the gas 
that is currently forming stars highlighted in green (see colorbar for relative values). The third and fourth rows from the top show the old stellar 
component in red, and newly formed star particles in black. The bottom panels show the gas contours with the gas velocity perpendicular to the 
contour plane. 
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Fig. A.2. Same as Fig. lA. ll for model A with Rg = IRj, with the c = 25 satellite on a planar, less radial orbit. 
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Fig. A.3. Star formation rates for model A with Rg = Rj during the 1:5 
merger with a satellite on a prograde, very radial orbit with an inclina¬ 
tion with respect to the plane of the disk of 0 or 30 degrees. 



time [Gyr] 


Fig. A.4. Star formation rates for model Cl with Rg = R^ during the 
1:5 merger with a satellite on a prograde, very radial orbit with an incli¬ 
nation with respect to the plane of the disk of 0, 30, or 60 degrees. 
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